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O Abstract 

In the Bondi formulation of the axisymmetric vacuum Einstein equations, we argue that the 
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"surface area" coordinate condition determining the "radial" coordinate can be considered as part 
^ — , . of the initial data and should be chosen in a way that gives information about the physical problem 

whose solution is sought. For the two-body problem, we choose this coordinate by imposing a 

\D 

■ condition that allows it to be interpreted, near infinity, as the (inverse of the) Newtonian potential. 
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In this way, two quantities that specify the problem - the separation of the two particles and 
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their mass ratio - enter the equations from the very beginning. The asymptotic solution (near 
infinity) is obtained and a natural identification of the Bondi "news function" in terms of the 



source parameters is suggested, leading to an expression for the radiated energy that differs from the 
5— i , standard quadrupole formula but agrees with recent non-linear calculations. When the free function 

bjr>. 

of time describing the separation of the two particles is chosen so as to make the new expression 



^ ■ agree with the classical result, closed- form analytic expressions are obtained, the resulting metric 

H ; 

approaching the Schwarzschild solution with time. As all physical quantities are defined with 
respect to the flat metric at infinity, the physical interpretation of this solution depends strongly 
on how these definitions are extended to the near-zone and, in particular, how the "time" function 
in the near-zone is related to Bondi's null coordinate. 
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I. INTRODUCTION 



The characteristic initial value formulation Jl] of the axisymmetric vacuum Einstein equa- 



tions, first proposed by Bondi and co-workers 



has two distinct advantages over the initial 



value formulation in terms of a space-like hypersurface jjj: (i) there are no constraint equa- 
tions to be satisfied by the initial data and (ii) the unknown functions can be determined, 
one at a time, by performing an hierarchical series of quadratures over the initial data and 
"known" quantities (including functions of integration from previous quadratures) j^, U, 0] . 
Thus there are no partial differential equations in three independent variables to be solved! 

This is achieved by writing the equations in terms of coordinates adapted to a series of 
null hypersurfaces u = const, one of which is the initial hypersurface (assuming that the 
spacetime, outside a world tube containing the source, can be foliated by such hypersurfaces), 
and requiring that all but one coordinate be constant on the null geodesies (generators) in 
these null hypersurfaces. Then the "main" equations are ordinary differential equations 
that propagate the initial data along these geodesies ("rays" or bicharacteristics). The 
radial coordinate - the one that varies along the rays - must be specified by imposing one 
additional condition. 

There are two standard choices for this additional condition, one of which is always 
made, resulting in the so-called "area" or "affine" radial coordinate, respectively. After 
either choice is made and the equations integrated, the dynamics of the two-body problem 
must be deduced by interpreting the properties of the two bodies (masses, separation) in 
terms of the functions of integration in the chosen coordinates. This is a non-trivial matter 
involving a lot of approximations j^, 0| . 

Now, for a unique solution, additional (boundary) data must be given on a timelike (a 
world tube) or null hypersurface intersecting the initial u = const hypersurface. Such a 
timelike boundary hypersurface is used in the "Cauchy Characteristic Matching" (CCM) 
approach 5.7, and is conventionally defined as having constant "radial coordinate", while 
in the interior (Cauchy) problem it is identified with a coordinate sphere (in a different 
coordinate system). 

In this paper we propose that the choice of radial coordinate should be made in a way 
that is appropriate to the geometry of the physical problem, i.e., that the surfaces of con- 
stant radial coordinate must be identified with the level surfaces of a physical quantity. 
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In this way, the definition of the radial coordinate implicitly defines the physical problem. 
Applying this idea to the two-body problem, we impose a condition that follows from the re- 
quirement that the two-surfaces of constant radial coordinate in any null hypersurface share 
a geometric property with the surfaces of constant Newtonian gravitational potential of two 
point particles in Euclidean 3-space. This choice has the advantage that important physical 
parameters, such as the separation 1 of the two bodies and their mass ratio, appear in the 
equations from the very beginning. They thus act as "source" terms for this problem. We 
then integrate the equations, assuming a series expansion near infinity and asymptotically 
flat boundary conditions. By considering the angular dependence of the "mass" parameter 
of the solution, the arbitrary function of two variables generating the solution is then nat- 
urally identified with the first time derivative of the quadrupole moment per unit mass of 
the source. 

Section II introduces the notation we will be using and section III describes the choice 
of radial coordinate. In section IV we integrate the Einstein equations assuming a series 
expansion of the unknowns near infinity. In section V we use physical arguments to determine 
the arbitrary function in the solution, arriving at a result which differs from Bondi's. We 
argue that this discrepancy is due to our asymptotic definitions of physical quantities. In 
section VI we discuss the validity of our results, pointing out that a similar choice for the 
radial coordinate might be convenient for investigating the near field also. Finally, in section 
VII, we state explicitly all assumptions in our approach, compare our derivation to Bondi's 
and clarify the conditions under which our calculations can be interpreted as modelling the 
real physical problem. 

II. FORMULATION OF THE PROBLEM - NOTATION 

We will restrict our considerations to the the case where the axial Killing vector is hy- 
persurface orthogonal (no particle spin). We will use the symbols u, £, i], ip for the four 
coordinates. The assumption that the Killing vector trajectories (parametrized by ip) are 
hypersurface orthogonal implies that g utp = g^^ = g 1]Lf> = 0, while the requirements that 
u is a null coordinate and the curves along which only £ varies are null geodesies imply 

1 The word "separation" does not mean "physical distance" . For the sense in which physical terms are used 
in this paper, see the Remark at the end of section III. 
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Q] that g uu = g ur> = also. Thus there are 5 remaining metric functions of u, £, rj, de- 
noted by the letters B, V, K, R, U. We will write the line element in terms of a canonical 
Newman-Penrose null frame 2 as ds 2 = I ® n + n ® Z — m (g> m — m £g> m, where: 

1 V -1 /K 2 \ 

' = ~B dU ' n = d£ + —du, m= -j=\—(dr]-Udu)+iRd^y 

a » f 9 V 9 tt 9 \ n 5 t 1 /ii a i 9\ 

The frame transformations preserving the Z direction and keeping m-m parallel to the 
Killing vector are a boost in the l-n plane and a null rotation around I parametrized by S 
and T, respectively: 

l^l/S, S[n + T(m + m)+T 2 l], m^m + Tl. (2) 



If one chooses S = H^, T = —RE )V / [y2K 2 H^) and redefines 5, V and C/ appropriately, the 
new frame is the canonical frame that is obtained by a redefinition of the radial coordinate 
£ — > S(w, £, 77). Thus an additional condition must be imposed to eliminate this freedom, 
effectively reducing the number of unknown functions of three variables from 5 to 4. Bondi's 
standard notation (and choice of radial coordinate, K = £) is obtained by making the 
substitutions: £ — > r, 77 — > — cos#, £> — > exp(— 2/3), V — > Vs/r, K —> r, U — > sin^f/^, i? — > 
r sin (9 exp(— 7). 

The reason for using £, 77 instead of the more usual r, x(= cos 9) is that we want to reserve 
the latter symbols for the radial and angular coordinates of a Euclidean prolate spheroidal 
coordinate system, which we will use to impose a condition on the function K(u, £, rj) in 
the next section. Here we note that, in such a coordinate system (with singularities at ±0 
on the z-axis), the Schwarzschild solution for a mass M at z = a on the symmetry axis is 
given by the substitutions: 

r x — a 2M 
4^r — ax, r) , B — > 1, V — > 1 



r — ax r — ax 

2 2\ /7i „2> 



/Y-^r-ax, Z7 — »• 0, i2 -> ^(r 2 - a 2 ) ^(1 - x 2 ). (3) 

In the limit a —>■ 0, we obtain the standard form in spherical r, 2 coordinates. What is 
more, J3J) remains a solution if we allow a to be a function, a(u), as it is simply a change of 



2 We use the standard NP notation the complex null-tetrad basis {I, n, m, m} is normalized to 
I n = —m • to = 1. We will restrict the symbols {I, n, m, rn} to denote the basis co- vectors (differential 
forms) while the vectors (differential operators) will be denoted by the standard symbols {D, A, 8, 5}. 



coordinates applied to the Schwarzschild solution in "spherical" £, r] coordinates: 

B = l, V = l-?j-, K = £, U = 0, R=tJ{l-r ] 2 ). (4) 

The 2-surfaces of constant u, £ plavan important role in the characteristic initial value 
formulation of the Einstein equations |l|,|5|,lZ|- I n particular, the induced metric on them, 

do* + #6^, (5) 

is known as part of the initial data. Thus imposing conditions on this part of the metric, in 
a sense, defines the physical problem whose evolution is determined by Einstein's equations. 
Of course, only the rj dependence of K is determined by the initial data, and this can be 
eliminated by a redefinition of the rj coordinate. This is the approach usually taken in 
imposing the "spherical area" radial coordinate condition K = £. We will instead impose 
a different condition that will allow these 2-surfaces to have different shapes depending on 
the ratio a(u)/£, where a(u) is a measure of the separation of the two particles. 

III. THE GEOMETRY OF THE SURFACES OF CONSTANT NEWTONIAN PO- 
TENTIAL 

Prolate spheroidal coordinates (in Euclidean 3-space) is the appropriate coordinate sys- 
tem to use for the two particle problem. The Euclidean line element in such a coordinate 
system with singularities at ±a on the z-axis is 

ds% LAT = (r 2 - aV) (-£L + + (r 2 - a 2 )(l - x 2 ) d^ 2 , (6) 

and the solution of the Laplace equation for the Newtonian potential of two point particles 
with masses mi at z = a and m 2 at z — —a is 

V N = - -!?-. (7) 

r — ax r + ax 

For the 1-particle problem (say 1712 = 0), Bondi's geometrical coordinate choice K — £ = 
r — ax, admits several interpretations in terms of properties of the surfaces of constant 
Newtonian potential: 

• K is constant on the surfaces of constant Newtonian potential 
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• 2/K 2 is the Gaussian curvature of the surfaces of constant Newtonian potential 

• 1/K A is the flat-space norm of the gradient (?7 a6 Viv,aViv,b) of the Newtonian potential 
of a unit mass. 

It is tempting to impose a condition, following from such an interpretation of K, on the 
metric functions in (JTJ), with (V — 1)/2 playing the role of the Newtonian potential. However, 
apart from ambiguities in relating the metric function V to the "Newtonian potential", 
involving the function B, one must remember that the function V is not a part of the data 
on the initial u hypersurface and, therefore, relating it to K for all £ (and u) may unduly 
restrict the evolution. We will instead impose a condition that will determine K(u, £, rj) 
without reference to the other metric functions. 

To do this we write the flat-space metric (jUJ) in terms of an orthogonal system of coor- 
dinates £, r), (p, where the £ coordinate is the (inverse of the) Newtonian potential (0) per 
unit mass: 

mi + m 2 mi m 2 , x rx — a rx + a 

= 1 , {m 1 +m 2 )i] = m 1 h m 2 . (8) 

4 r — ax r + ax r — ax r + ax 

Let us denote by iV2(0, ip) the inner product of the gradients of the functions 0, ib with 

respect to the flat metric 

»r„/ , ,\ (r 2 — a 2 ) <t> r ib r + (1 — x 2 ) <t> x ib x 

iV2(0, = — ; 9 V . — — — • (9) 

r z — a x 

Then it is easy to verify that £, 77 as functions of r, x given by (JSJ) satisfy A^2(£, 77) = 0. It is 

remarkable that orthogonality is maintained when the "radial" and "angular" coordinates of 

the individual 1-particle problems [see @] are combined linearly in this simple way. This is 

undoubtedly due to the appropriateness of prolate spheroidal coordinates for describing this 

problem. In terms of a parameter A depending on the ratio of the masses, the coordinate 

transformation (jHJ can be written 

r 2 — a 2 x 2 a (Ar + ax)(l — x 2 ) rti\ — m 2 , . 

i = — — , 7] = x — , where A = ■ . 10 

r + Aax r z — a z x z m\ + m 2 

We note that, when x — > ±1, also rj — > ±1, while as r — > 00, £ — > r and 77 — > x. Thus, for 

r ^> a, the coordinates £, rj behave as prolate spheroidal (or spherical) coordinates r, x. 

In terms of the new coordinates defined in the flat metric (JHJ) takes the form 



Now, it is easy to verify that N2(i], 77) = iV2(£, £)(r 2 — a 2 )(l — x 2 )/£ 4 , so that the metric 
on the £ = const surfaces in Euclidean 3-space can be put in the form (jSJ) with 

R 2 = (r 2 - « 2 )(1 - x 2 ), ^ 4 = = where (12) 

a 2 (l - A 2 )(-r 2 + 3r 2 x 2 + 4 Xarx 3 + a 2 x 2 + X 2 a 2 x 4 ) 
S= {r + XaxY ' (13) 

The denominator of S (or £) cannot vanish as the coordinate ranges are r > a, —1 < x < 1, 
and the parameter A is by definition less than one in absolute value. 

We are now ready to make our choice for the Bondi "radial" coordinate £ that is appro- 
priate for the two particle problem: we will require that K will be given by (|12j) for any 
value of u, i.e., that the parameter a in (fT3*j) . determining the separation of the two particles, 
will be allowed to be a function of u. 

To be used in Einstein's equations K must, of course, be given as a function of u, £, 77 
by inverting the coordinate transformation (jlOjl . This cannot be done in closed form, as it 
depends on the roots of a 5th order polynomial in £, 77. However, it is easy to convert the 
series expansion of K for r — > 00 to one for £ — > 00, making use of the limits noted after 
(110}) . This is most easily done by replacing, in the series expansion of S with respect to r, 
the symbols r, x by £, 77, respectively, one order at a time and expanding again: 

s _ (1 -X 2 )(-l + 3x 2 )a(u) 2 _ AX(l-X 2 )x(-l + 2x 2 )a(uf + 1 4 



(l-A 2 )(-l + 3r/ 2 )q(w) 2 _ 4 A (1 - A 2 ) x (-3 + 5x 2 ) a( M ) 3 + ^l^ 4 
(l-A 2 )(-l + 3r7 2 )a( M ) 2 4 A (1 - A 2 ) 77 (-3 + 5T7 2 ) a( M ) 3 4 4 

« + « ° H • 



Substituting in (fT2"|) we obtain, to this order, 



K „ ^ _ (l-A 2 )(-l + 3r / 2 )a( M ) 2 + A (1 - A 2 ) rj (-3 + 5 V 2 ) a(u) 3 + Q{ 1 } 3 _ ^ 



In this way we can compute the series expansion of S, and therefore of K, as a function 
of u, £, 77 to any desired order. We note that the angular dependence of the first two 
"correction terms" in K are proportional to the Legendre polynomials P2, -P3 and therefore 
vanish when integrated over angles. Thus we can interpret our coordinate £, alternatively, 
as an approximate "luminosity distance": the area of the surfaces of constant u and £ equals 
4 7r £ 2 to a good approximation. To find out how good this approximation is we must carry 



the expansion to 5 more orders! We then find (see (|5|)): 



JdS = j K^dp = 4*? (l ~ |(1 - A 2 ) 2 (l - 6 A 2 )^ + 0{± 



7\ 



(16) 



'£ J 

Thus our radial coordinate £ is also a "luminosity distance" to a very good approximation. 

It should be pointed out that the reflection symmetry, 77 — > —rj, assumed in many inves- 
tigations 0,0], holds here only if, at the same time, A — > —A. This is to be expected, as the 
interchange of the positive and negative z-directions will not give the same problem unless 
the two particles are also interchanged. 

Finally, we compute the flat-space quadrupole moment of the two particles relative to 
their center of mass. In terms of cylindrical coordinates p, z, p with origin at the center of 
mass, the mass density distribution is 

p = ^A{ mi S[ z - (1 - A) a(u)\ + m 2 5[z + (1 + A) a(u)]}. (17) 

2 7T p 

Then the quadrupole moment tensor can be written 

D ab = J p{3x a x b -x 2 5 ab )dV = Q{u){3 5 a z 5 b z -5 ab ), (18) 

where 

Q(u) = J {m^lz - (1 - A) a(u)] + m 2 5[z + (1 + A) a(u)]} z 2 dz 

= (mi + m 2 )(l-A 2 )a(n) 2 , (19) 

and we have used the definition of A to write m 1? m 2 in terms of mi + m 2 and A. 

Remark: All statements regarding physical quantities (particle separation, mass density, 
quadrupole moment, etc.) in this paper refer to an unphysical background Minkowski space, 
which is the limit as £ — ► 00 of the sought solution of Einstein's equations. In the physical, 
curved, spacetime they are to be understood only as names referring to the combination of 
variables in their definition. Thus, for example, "particle separation" means 2a(u) and not 
physical distance between the two particles. 

IV. SERIES SOLUTION OF THE VACUUM EINSTEIN EQUATIONS 

Knowing K and assuming a formal series expansion for the function R, valid near £ — > 00, 
where R approaches the flat-space limit (@J), we can carry out the well-known hierarchical 
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series of £ integrations to obtain the other metric functions. Thus, we assume that as £ — > oo, 
R has the series expansion 

^(^-fl(l + ^^+^+0(|)'). (20) 

Then the Ricci component $00 i n the frame defined in (0) can be solved for B^/B, giving 

B 1 = _K 1 + 2 R 1 _K I1 _KRl 

B K R R 2 K^ K ' 

Integrating the series expansion of the rhs when K, R are given by equations (|15j) . fl2()L and 
choosing the function of u, 77 of integration to satisfy the boundary condition (0J at infinity, 
we find 3 



b ~ 1 + v ' M ; ^ 1V "' ^ + . . . (22) 



(l-377 2 )(l-A^a(«r + 2 Cl (u, 77) 
4£ 

We next compute the series expansion of the Weyl tensor component ^0, obtaining 4 



(1 - 3fy 2 ) (1 - A 2 ) a(u) 2 + 2 Cl (n, r?) 2 -4c 2 («, 77) 

W - r • • . (,A5J 

Now, the condition for asymptotic flatness (absence of incoming radiation at infinity) is that 
4I ^0 — 0(|) 5 , so that we must require that 

C2( „,,) = ( 1 - 3 " 2 )( 1 - A2 );W 2 + 2 ^") 2 . ,24) 

Next, the Ricci component $01 can be written 

[EJ^l] = K 2 [^-^ + 2^ + 2^]+2i^[^-2^], (25) 
[ R 2 \ [ B 2 B R 2 R i[ B R h K J 

so that, with the series expansions of K, R, B known, U can be obtained. We find 

[1 - T] 2 ) Ci tV (u, 7]) -2r]d(u, 77) 



e 

[1 - r] 2 )[A ci(u, rj) c hv (u, 77) - x (u, 77)] - 8 77 c x (u, rj) 2 + g 



3£ s 



3 We only give the leading two terms in each equation. All calculations were carried out using Mathe- 
matica and the author's "RGTC" package (http://www.inp.demokritos.gr/~sbonano/RGTC/), keep- 
ing terms to one order higher than in JT5J, (|2U|) . i.e., to four orders beyond flat space. The 
complete expressions can be found in the Mathematica notebook "NewRadialCoord.nb" , available at 
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http : //www. inp . demokritos . gr/~sbonano/TwoBody/NewRadialCoord.nb 

In the frame defined in jy, all NP quantities (spin-coefficients, Ricci and Weyl tensor components) are 
real - for a proof see the Appendix in [10(. 



9 



where xq(u, rf) is an arbitrary function of integration. The second such function, to be 
added to U, must be set equal to zero for asymptotic flatness. 

Proceeding in the same way, equation $ n + 3A = can be solved for (V K K^) % in terms 
of known quantities, and, with V known, equation $ 2 — gives (K R u / R) ^. Denoting by 
yo(u, rj), z (u, rf) the functions of integration, V and R tU are given by 

V ~ 1 + + ... (27) 

R >u ~yf(l- rf) [z {u, rf + 2ci(u,^) 20 (u, q )+(i-3^)(i-A^) o (u)d( u ) j + _ _ _ (2g) 

Finally the requirement that this agrees with the u derivative of (}2T)j) determines the u 
derivatives of the coefficients rf) (except for 02(1*, rf) which is given by (J2H)): 

ci, tt (u, rf) = z {u, rf), c 3 , u (u, rf) = . . . (29) 

This completes the integration of the so-called "main" equations. Of the remaining equa- 
tions, $11 = is satisfied identically, while the vanishing of $12, $22 impose the following 
two conditions - conservation laws - on the three functions of integration xo(u, rj), yo(u, rf) 
and z (u, rf = c ljU (u, 77): 

Xo,u = -2/0,7? + cic hur] -3c 1)U Ci tV + 3 r] (1 - A 2 ) a(u) a(u) 

y 0>u = 2c\ u -2c 1)U + (1 - r] 2 ) c liUTIV -4:r]c ljUV (30) 
+ (1 -3t7 2 )(1 - \ 2 )\a(uf + a{u)a{u)\. 

The series solution obtained in this section reduces to the one with Bondi's coordinate 
choice K = £ in the limit A 2 — ^ 1. In fact it can be obtained as a coordinate transformation 
from the Bondi solution if one redefines the q to be given in terms of the corresponding q b 
(= coefficients in the series expansion of R in powers K) by the expressions obtained when 
K is replaced by its series expansion ([T5|). One then finds that C\ — C\b and c 2 is given by 
(}24|) . To make the conservation laws (|30j) take the Bondi form, one must also redefine the 
functions of integration xq(u, rf), yo(u, rf) as follows: 

x {u, rf) -> x 0B (u, rf) + | rj (1 - A 2 ) a(u) 2 , (31) 
y (u, rf) -> y QB (u, rf) + (1 - 3 r/ 2 )(l - A 2 ) a(u) a(«). 

Despite the fact that, formally, the solution obtained here is but a coordinate transformation 
of the Bondi solution, the explicit appearance of the "source" terms involving a(u) and A 
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in the conservation laws ([30)1 makes the choice of the arbitrary function c± describing the 
two-body problem and, consequently, the physical interpretation of the solution particularly 
simple and transparent. 



V. PARTICULAR SOLUTION DESCRIBING THE TWO-BODY PROBLEM 

The solution obtained in the previous section is a "general" solution of the equations in 
that it depends on the arbitrary function ci(u, rf). To obtain the solution for a particular 
problem, an appropriate choice for this function must be made. With the coordinates chosen 
to fit the two-body problem, it is reasonable to expect that the required arbitrary function 
will have a simple form. First, non-singular behavior on the axis (rj = ±1) requires that 
Ci(u, rf) = q(u, rf) (1 — rf) for some q(u, rf) that is well behaved at r\ — ±1. Making this 
substitution in the second conservation equation (|3T)J) and replacing yo(u, if) by — 2M(w, if) 
(Bondi's "mass aspect" definition - see equation (|27jl ). we obtain 

2M, a = -2(\-rff q, a 2 
-(1 - 37? 2 ){(1 - A 2 ) \a(uf + a(u) a(u)} - 4q, u } (32) 

-(1 - i] 2 )[(i - i] 2 )q,u VV -8vq,ur,}- 

Now, the rhs of this equation describes the energy loss of the system. On physical grounds, 
it must be negative definite and have the angular dependence that is appropriate to the 
one-dimensional motion of the two particles. This angular dependence, being independent 
of £, can be identified with the angular distribution of the flow of energy at infinity obtained 



in the linearized theory using the Landau-Lifshitz pseudotensor (see ll| equation 110.15) 



„ {±(D ab n a n b ) 2 + l -D ah D ah - D\n b D ac ii% (33) 

where D a ^ is the third time derivative of the quadrupole moment tensor 5 (J18)) and n a are 
the Cartesian components of the unit vector in the direction of propagation, which in our 



5 It is consistent to use the flat-space definition i|18|l for the quadrupole moment tensor here, as the L-L 
pseudotensor is defined in terms of the linear solution to the field equations, which is determined by a 
flat-space integral over the sources. Besides, only the angular dependence of D a b is used in the present 
argument, not the exact form of Q(u). 
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coordinates equal (J (1 — rj 2 ) cos<^, y (1 — ?7 2 ) siny9, 77) at infinity. We find 

^[1(3^-1)' + 1(6) -(3„°+l)] 

= l^^f + l-(3,f-l)] = l( l f-lf- (34) 

Thus the angular dependence of the mass-loss equation (}3*2*)) will be ~ (1 — r] 2 ) 2 , as appro- 
priate for this system, if we choose 

q (u, 7]) = |(1 — A 2 ) a{u) a(u), so that (35) 
d(«, 17) = |(1 - A 2 ) a(u) o(«) (1 - r? 2 ) = §Q (1 - r? 2 ), (36) 

where Q = (1 — A 2 ) a(u) 2 is the quadrupole moment per unit mass - see (|19jl . The same 
result follows from the requirement that the rhs of ()32|) be negative for all u, r], so that 
energy flows out of the system in all directions and at all times: (j33j) is the unique solution 
that is regular on the axis and makes the linear (and "source" terms) on the rhs of ()32|) 
vanish. 

With this c\(u, rj), the first conservation equation can be integrated giving 

x (u, V )=V(1- A 2 ) a{uf [\ + (1 ~ 7/2)(1 ~ A2) ^ )2 ] + 2/ M, v du, (37) 
where M(u, rj) is determined by the equation 

= -d - , 2 ) 2 ( (i - A2)KM) ; +a(M)a(M)1 ) 2 - -d - (i) 2 (as, 

once the function a(w) is known. 

Finally, using Q, the quadrupole moment per unit mass of the system, the leading terms 
in the components of the Weyl tensor in the frame defined in with ci(u, rj) given by 
(JHBJ), are: 

*o = |r h(u, V ) + Xr ] (3-5r ] 2 ) Qa{u)\ + . . . , 
^-f(M+^f^ + ^M)+..., (39) 

Using the evolution equations (|2*9|l . it can be shown that \I/o is proportional to (1 — rj 2 ). 
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We note that equation (j3o|) . giving the radiated energy directly in terms of the parameters 
A and a(u) describing the source, is an exact result following from a particular, physically 
motivated, choice of the arbitrary function Ci(u, if) (or q(u, if)). And this choice of q(u, rj), 
which eliminates all linear terms on the rhs of the energy-conservation equation (|32j) . is made 
possible by the existence of the extra term introduced by our choice of the function K, which 
fortuitously has the correct angular dependence (1 — 3i] 2 ). The final expression for the mass 
loss (J58|L relating the radiated energy to the square of the second "time" derivative of the 
"quadrupole moment per unit mass" differs from the standard quadrupole formula (third 
time derivative of the quadrupole moment) obtained in the linearized theory. However, 
it closely resembles the results of approximate non-linear calculations, where the radiated 
energy is found to be proportional to the square of an integral with respect to u of Q/M 
weighted by exp (2u/M) HQ, where M is a large constant "background" mass. Thus 
our choice of ct(w, rj), equation (f3*fi|). implying that M tU ~ Q 2 in the fully non-linear case, is 
in remarkably good agreement with these results. The discrepancy with the classical result 
must be attributed to the different definitions of "time" and "quadrupole moment per unit 
mass" (see the Remark at the end of section III). And the following considerations reinforce 
this conclusion. 

Bondi's approximate linear calculation for the news function, which agrees with the con- 



clusions of linearized theory 11 1 and post-Newtonian calculations lj,ll5|, gives Ci(u, rf) 
— ^Q(u) (1 — if) (in our notation Bondi's c equals —cf). If we want to reconcile the two 
results we must require that, as a consequence of the equations of motion 6 , the quadrupole 
moment satisfies the equation Q{u) = —\Q{u)/{m\ + mf) so that, 

Q(u) = Q exp[ — — -], and therefore, 

4 (mi + m 2 ) 

a{u) = a exp[-— - — -^]. (40) 
8 (mi + m 2 ) 

Of course, a(u) should be determined by the equations of motion following from the vanishing 
of the divergence of the energy momentum tensor of the two particles. But this requires 
knowledge of the field in the vicinity of the particles, which is well beyond the scope of the 
approximate calculations near null infinity carried out here. 



6 The equations of motion are used in the derivation of the result in linearized or post-Newtonian theory 



t rie eqr 



and, implicitly, in Bondi's derivation. 
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Nevertheless, we point out that, with a(u) given by ()4Uj) . all remaining u integrations (see 
equations (j2HJ), (j3*Tj) . (J3~%)0 can be evaluated analytically in terms of elementary functions 
(exponentials); and choosing the functions (of if) of integration to vanish (except for the 
final mass), the resulting solution will smoothly approach the Schwarzschild solution as 
a(u) — > exponentially with u — > oo. Now, if we assume that the Minkowskian definitions 
of "separation" , "velocity" and "time" used here approximate adequately the corresponding 
physical quantities, this solution has the unphysical feature that the two particles approach 
each other at a diminishing rate (both a(u) and a(u) decrease with time): the collision seems 
to end with hardly a whimper rather than a bang! However, the physical interpretation of 
"velocity" depends on the definition of "time" : near the source the time function, which can 
be taken to equal u + £ at infinity, is expected to have a logarithmic singularity, as in the 
Schwarzschild case 16|, 17]. A time function of the form t = u + ^ — (1 + k)a(u) + 8 (mi + 
m 2 ) log(£/a ) (where fc=const) evaluated on the curve £ = a{u) gives t = u — ka(u), so that 
da/dt = —1/k. Thus, with a proper definition of time in the near-zone, the rate of approach 
implied by ()4()|) might be physically acceptable. It seems too speculative to discuss further 
this possibility. But this example shows that one cannot simply identify w-derivatives at 
infinity with time-derivatives in the near-zone, and that the mass-loss formula ()38|) may be 
consistent with the standard quadrupole formula! 

In any case, this particular solution, being an exact analytic solution (in the form of a 
series expansion which can be continued to any order) of the Einstein equations, including 
the evolution equations, which approaches the Schwarzschild solution as ti — > oo, can be 
useful in testing the accuracy of numerical codes. For this reason, explicit expressions for 
the metric functions and the Weyl tensor components, together with the verification that 
the metric is Ricci-flat to the appropriate order and becomes Schwarzschild as u — > oo, are 
included in the Mathematica notebook "NewRadialCoord.nb" mentioned in footnote 3. 



VI. DISCUSSION 

The freedom in the choice of coordinates inherent in general relativity is invariably used 
to simplify the equations. This is due to the complexity of the equations. However, a 
coordinate choice that best simplifies the equations may not describe the physical problem 
in the most natural way. For example, in the static, axially symmetric problem, Weyl's 
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canonical p, z coordinates require that a physically spherical source be described as a linear 
distribution of mass. 

In the Bondi formulation of the axisymmetric vacuum equations, the essential simplifi- 
cation comes from the use of a null coordinate u and angular coordinates rj, tp which are 
constant on the null rays. The choice of parametrization on these rays does not simplify the 
mathematical problem any further. In this paper we propose that it should be made in a 
way that gives information about the particular physical problem whose solution is sought. 
For the asymptotic solution near infinity of the two-body problem, imposing a condition 
that follows from interpreting the surfaces on which u, £ are constant as surfaces of constant 
Newtonian potential, allows us to relate the arbitrary function of two variables generating 
the solution to two parameters describing the source, A and a(u), by a symmetry argument 
without making any approximations. This should be compared to Bondi's approximate and, 
in his own words, "distinctly crude" 7 derivation of the corresponding result. 

Despite its simplicity and directness, our derivation of the form of c\ is subject to a 
serious objection, as the result seems to depend on a choice of coordinates. How can such 
an arbitrary choice lead to a physically meaningful result? There are two answers to this 
criticism: (i) Once the arbitrary function in the solution (q(u, n)) has been related to the 
arbitrary function (Q(u)) determining the coordinate transformation (|15J). the latter is no 
longer arbitrary. The same result could have been obtained with Bondi's choice of radial 
coordinate = K, had one chosen the arbitrary functions xo(u, 77), yo{u, rj) to include extra 
source terms as in equation (jHljl . and then transformed the radial coordinate from K to £ to 
eliminate the extra term in the "mass aspect" 8 (this is discussed further in the next section). 
In this sense, the leading term in the coordinate transformation (jl5j) is determined by the 
solution, (ii) Our choice of radial gauge, following from the identification, near infinity, of 
the two-surfaces of constant u, £ with the surfaces of constant Newtonian potential, is not an 
arbitrary choice of coordinates but is closely related to the physics: the surfaces of constant 
£ are the wave fronts of the radiation emitted from the system, carrying information about 

7 See comments after equation (91) in An outline of Bondi's derivation is given in section VII. 

8 This extra term in the "mass aspect" does not change the "Bondi Mass", as its integral over angles 
vanishes. In fact, in obtaining his formula for cib, Bondi also obtains a first-order expression for Mb that 
has precisely the form f(u)(l — 3 cos 2 9), as in (fSTft : in going from eq. 88 to eq. 89 of reference Q, Bondi 
has set his function p(u) equal to 2Q o (plus a constant MO)! 
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its properties; and identifying their constant time sections (t = u + £ near infinity) with the 
surfaces of constant Newtonian potential gives a description that matches the dynamics more 
accurately than the use of either spherical or prolate-spheroidal [6] surfaces (any change in 
the dynamics at infinity must be manifested through changes in the Newtonian potential). 
In fact, the dynamics can be better approximated if one allows the parameters m\, m 2 to 
depend on u to reflect the relativistic velocity dependence of mass. All £ integrations remain 
unchanged and only equations involving u derivatives of the unknown functions will acquire 
extra terms, leading to more complicated evolution equations (|29|) and conservation laws 
(jHUj) . Indeed, the energy loss equation (jSHJ), expressed in terms of Q, remains unchanged 
even when the parameter A in Q is allowed to depend on u. 

This leads to the suggestion that, in the inner problem also, a coordinate condition 
based on a property of appropriately defined constant-Newtonian-potential (or constant- 
Newtonian-energy) surfaces should be used, as these surfaces have the topology of the "pair 
of pants" picture of the horizon near the source, the set of points with constant Newto- 
nian potential consists of two disconnected subsets, one around each particle. This can best 
be seen if, near £ = 0, one switches to the r, x coordinates used in the definition of £, rj fllOl) . 
Of course, near the source the null coordinate u will not be well behaved and a different 
time coordinate must be used, while care must be exercised in defining the "Newtonian 
potential" . 

VII. SUMMARY AND CONCLUSIONS 

This paper addresses the mathematical problem of solving Einstein's equations with the 
symmetries (one hypersurface-orthogonal Killing vector with closed spacelike orbits) and 
boundary conditions (asymptotically flat, no incoming radiation) that are appropriate for 
the two-particle collision problem. All calculations are carried out near infinity, and physical 
quantities are defined with respect to the limiting Minkowski space (see the Remark at the 
end of section III). The Bondi formulation of this problem is used as it allows free initial 
data and a sequential integration of the equations. Both the solution and the coordinate 
system are only assumed to be valid outside some sufficiently large time-like world tube. 

A particular, explicit, solution to this mathematical problem depends on two arbitrary 
choices: (a) the "radial" coordinate condition - the coordinate expression for the determinant 
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of the metric on the two-surfaces of constant u, £ (the function K(u, £, if)), (b) the free 
function c±(u, rj). Mathematically, these choices are completely arbitrary (modulo boundary 
conditions) . 

The basic idea of this paper is to try to "guess" the appropriate form for these arbitrary 
functions based on the "expected" properties of the solution describing the two-particles; 
properties which can be calculated in the region where the solution is valid, i.e., outside 
the large world tube. Thus, (a) the function K(u, £, rj) is chosen by the requirement that, 
on any surface of constant u and £, it is identical to the corresponding function in the 
metric of the surfaces of constant Newtonian potential of two point particles in Euclidean 
3-space. Even though the expression for K in Euclidean space is valid everywhere, we 
only use the asymptotic form of K, valid for, say, £ > 10a(w). (b) The function Cx(u, rj) 
(or q(u, 77)) is chosen by the requirement that, at infinity, the angular dependence of the 
emitted radiation is that appropriate to quadrupole radiation. The same q(u, rj) satisfies 
the physical requirement that a non-negative amount of energy is radiated in any direction 
and at any time. 

The physical significance of allowing K to have a term f(u, rj) /£ is that it changes the 
"mass- aspect". This is most easily seen by making the change of coordinates r — > r + 
p(u, 9)/f to the Schwarzschild metric in null-spherical coordinates, ds 2 = (1 — 2M/r)du 2 + 
2 du dr—r 2 dQ 2 . Then, to first order in 1/f, the metric acquires an extra term 2dud9pe/f but 
the "mass aspect" (=coefficient of 1/f in g uu ) is also changed to M—p u due to a contribution 
from the 2 du dr term in the metric. This is the basic reason for the existence of the extra 
term on the rhs of the mass-loss equation for the metric defined by (jljl . which allows the 
linear terms to vanish, leaving M = — c'i 2 , when one chooses K = ^—f(u)(3r] 2 — l)/^+. . . and 
C\ = — rj 2 ). One could then make the ad-hoc choice f(u) = Q(u) and claim that one 

has obtained the standard quadrupole formula as an exact result of the non-linear equations! 
But this ad-hoc choice lacks any physical justification for calling this Q the "quadrupole 
moment" . Our interpretation of £ in terms of the Newtonian potential at infinity gives both 
the correct angular dependence to the l/£ term in K and provides a physical interpretation of 
the coefficient f(u), namely Q(u) / (4M) (with Q, admittedly, defined in a flat background). 
Of course, the question of uniqueness remains: could a different physical interpretation of 
the "radial coordinate" £ lead to a function K that (i) agrees with Bondi's K — £ both at 
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infinity and for the 1-particle limits (A — > ±1, or a(u) — > 0), (ii) the limit lim £(K — £) 
exists and has the form f(u)(3n 2 — 1) but with a different physical meaning of f(u)7 

It is appropriate, at this point, to recall how Bondi obtains his relation between c\ and 
the quadrupole moment. He begins with the static, axially symmetric vacuum metric in 
Weyl coordinates 

in which the function if> satisfies the flat-space Laplace equation. He then 

1. Takes for ip the solution (in spherical coordinates p = .RsinG, z = -RcosO) 

_m_L>cos6 (Q + l/3m 3 )(3cos 2 - 1) 
^~~R R? ' 

m, D and Q being called the "mass", "dipole", "quadrupole moment" of the source. 
The term 1/3 m 3 is needed to make Q = for the Schwarzschild solution, as it is 
obtained from a linear mass distribution in these coordinates. 

2. Finds the coordinate transformation to Bondi coordinates (t, R, 6) — > (u, r, 6) as a 
series approximation for large r (or R). 

3. By comparing coefficients of the transformed (to Bondi coordinates) Weyl static solu- 
tion with the series expansion of the metric (0), he finds expressions for the unknown 
functions C3, Xq, (equivalent to his C, N, M) in terms of m, D, Q. The function c\ 
is also obtained as a coordinate-dependent expression. 

4. Next he supposes that, for slowly varying fields, he can use the same expressions for 
these functions allowing the constants m, D, Q to be functions of u. And, for weak 
fields, he solves the linearized form of the equations (|29jl. (f3T)J) determining the u 
derivative of C3, Xo, yo (expressed now in terms of m, D, Q) for the independence of 
D, m, ci, respectively. In this way Q, instead of c\, becomes the free function. 

It will be observed that the last step involves two serious approximations. Bondi is 
alluding to this in calling his derivation "distinctly crude" . In particular, what conclusions 
would have he obtained had he allowed the constants m, D, Q to be time-dependent before 
making the coordinate transformation, and, perhaps, having found first a time- dependent 
first-order correction to the static solution? This crucial step (allowing the constants to 
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be functions of u) is necessary in order to generalize the static Newtonian potential to a 
dynamic one, but it is also the weakest point in the derivation: when ip is ^-dependent, the 
Weyl metric is not a solution (to first order in rh, D, Q) of the equations! 

In our approach, we make contact with a u-dependent Newtonian potential via our inter- 
pretation of £, which we are free to do as it is a coordinate condition, not an approximation 
for slowly varying fields. Moreover, we do not need to integrate the linear approximations of 
(J2HJ), (jHUj) to relate the free function c\ to the properties of the source. We obtain c\ directly 
by a physical argument: it is the unique, non-singular solution that makes the linear terms 
on the rhs of (jH2J) vanish, resulting in an exact expression for the mass loss (M = -c'i 2 ) that 
has both the correct angular dependence and is negative definite. If we accept these argu- 
ments as physically reasonable, then we are obliged to conclude that the radiation depends 
on the "second derivative of the quadrupole moment per unit mass" - a result which seems 
to disagree with the standard quadrupole formula. We must remember, however, that our 
physical quantities (and those of references 12, are defined at infinity and the "dots" 



in our formula cannot be translated directly to "time derivatives" in the near-zone. On the 
other hand, the quadrupole formula has been proved only in a linear (in the metric 
tensor) approximation to the field equations. 

The entire solution depends on an arbitrary function of one variable, a(u), that is qual- 
itatively related to the particle separation. This function is determined by the equation of 
motion which, however, cannot be evaluated at infinity. 

In conclusion, we wish to stress that the solution obtained in this paper (with the choice 
(JHBJ) for ci) is a formal asymptotic approximation (for large £) to an exact solution of the 
full non-linear Einstein equations that has all the required properties for describing the two 
particle collision problem. The solution can be continued to any desired order in the expan- 
sion and appears to converge for £ > a(u) in the sense that no large numerical coefficients 
appear. However, it is safer to assume that it is only valid for, say, £ > 10 a(u). The inter- 
pretation of the solution, which was also the motivation for making the above choices for the 
free functions K and c\, avoids mathematically unjustified steps in Bondi's interpretation; 
however, it involves physical quantities which are defined only in a background Minkowski 
space - the limiting form of the solution as £ — > oo. Using this solution to calculate the 
outer boundary data for the inner problem in a CCM approach , the physical definition 
of the function a(u) and of "time" in the near-zone can, in principle, be obtained. 
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Note added in proof: It can be shown that the form of K given in ()15|) is valid for any 
mass distribution. Specifically, if we use spherical coordinates r, x and define £ (r, x) by the 
series expansion (Newtonian potential of a unit-total-mass distribution) 

11 Afi(u)ar M 2 (w)(3x 2 -1) M 3 (u) x (5 x 2 - 3) 

and r](r, x) by the requirement that V£ • V77 = 0, then, proceeding as in Section III, we find 
that 



( ^(^(3^-1) fc 2 ( M )7y(5 7y 2 -3) 
* = * + J + ^ + 



where k^u) = [M x {uf - Af 2 (u)]/4 and fc 2 (u) = -M 1 ( y u) 3 + lM 1 {u)M 2 {u) - |Af 3 (u). Thus, 
when the dipole-moment vanishes, fci(w) = — M 2 (m)/4 = —Q/(AM) for any mass distribu- 
tion. 
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